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We consider the dynamics of vortex strings and sound waves in superfluids 
in the phenomenological Landau-Ginzburg equation. We first derive the vortex 
equation where the velocity of a vortex is determined by the average fluid velocity 
and the density gradient near the vortex. We then derive the effective action for 
vortex strings and sound waves by the dual formulation. The effective action might 
be useful in calculating the emission rate of sound waves by moving vortex strings. 
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Recently there was some new interest in understanding the dynamics of vor- 
tex strings and sound waves in superfluids or superconductors [1,2]. While many 
aspects of the vortex dynamics has been analyzed extensively and used to many 
applications, we feel that some aspects of the basic vortex dynamics, like the vor- 
tex equation or the interaction between vortices and sound waves, are not well 
understood as one hoped. Based on recent progress in rigorous understanding of 
the vortex dynamics in relativistic field theoretic models [3], where the Magnus 
force on the vortex string is important, we derive hopefully a new vortex equation. 
While this is not a 'dynamical' equation of motion, it leads to some insight into the 
vortex dynamics. The model we use to describe superfluids and superconductors 
is the nonlinear Schrodinger equation studied before [4] . In addition, following the 
dual formulation of the theory of a complex scalar field [5], we derive an effective 
action for vortices interacting with sound waves. While many aspects of our goal 
have been studied before, there still seems to be a room for improvements from 
our point of view. 

In classical hydrodynamics a vortex moves with the local velocity of the fluid. 
Such a view has been applied to vortices in superfluids [1]. The quantization of 
the vortex strings has been also investigated [6]. However, the vortex equation we 
will derive has an additional term which depending on the spatial dependence of 
the fluid density and has an interesting physical interpretation. The vortex string 
equation we derive can be regarded as a nonrelativistic limit of the relativistic string 
equation. The relativistic vortex string equation has been useful in understanding 
the interaction between vortex strings very close to each other [3]. 

Since the work by Kalb and Ramond, there has been many works relating the 
antisymmetric tensor field coupled to vortex strings to vortex strings in superflu- 
ids and superconductors. However, the antisymmetric tensor field have usually 
described massless Goldstone bosons moving at speed of light rather than sound 
waves of speed less that unity [7]. The Kalb- Ramond antisymmetric tensor field 
is generated by the dual formulation of the massless Goldstone field. The nonzero 
fluid density appears as a uniform 'magnetic' field for vortex strings and give rise to 
the Magnus force on strings. Here we use the dual formulation to get the effective 
action for vortex strings and sound waves, extending the two dimensional result in 
the first paper of Ref. [5]. Recently, there has been also some interesting progress 
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in precise understanding of relativistic vortex dynamics where the Magnus force is 
important [8]. 

We start with the nonrelativistic Lagrangian for a complex scalar field, which 
may be coupled to the electromagnetic field, 

£ = - _ — |(V + ieA)^| 2 - ^(|^| 2 - p ) 2 + •••• (1) 

2 2t72- 2 

where the dots denote terms not important for our discussions like e^ol^l 2 . The 
average charge density is given by a constant po- The field equation is then 

i*= — — (V + ieA) 2 ^ + ... (2) 
2m 

There is a global Abelian symmetry whose charge density and current are p = \^\ 2 
and 
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^*(V + ieA)m - (V + ieA)^*^ 
The conservation law would be p + V ■ J = 0. 



(3) 



There are quantized vortex strings in the system. With ^ = ^fpe l6 ', the phase 
6 changes by 2tc around each vortex string. The natural core size of vortex strings 
would be given by the coherence length £ = 1/ \fgmpQ. To understand the vortex 
string dynamics, let us first consider a single vortex string whose position is given by 
R(t, a) with string coordinate a. The tangent vector along the string would be then 
R'(t, a) = dH/da. We choose the a sign so that the 6 field increases with the right- 
hand rule. We introduce two unit vectors Ni(t, a), N2(i, a) perpenducular to R' 
At the point on the vortex line, R', Ni, N2 would form a right-handed orthogonal 
basis for three space. There is a new coordinate system {a, C'X 2 ) 011 this basis 
near that point such that 

x = R(t, a{t, x)) + C A (^ 5)N A (t, a{t, x)) (4) 

^1=1,2 



Near the vortex, the relation between x and (a, C 1 , ( 2 ) would be invertable. 
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By evaluating the time derivative of Eq.(4), we get that on the point on the 
vortex 

C4(*,R) = -R-N4 (5) 

Similarly, the space derivative d/dx^ on Eq.(4) leads to 

R' 

W = 1W (6) 
V( A = N A 

on the point. Two spatial derivative V 2 of Eq.(4) leads to 

vV = (7) 

To obtain a vortex equation from Eq.(l), we first notice that near a vortex 
string, the scalar field would be ^ ~ tp s = (C 1 + i( 2 )- In general the scalar field 
can be written as \& = ip s &, where $ would not in general vanish on the vortex 
string. Let us evaluate Eq.(l) on a point on the string, to get 
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$V 2 ^ S + 2V$ • Vip s + 2ie$A • Vip s ) 



(8) 



The rest of terms vanishes on the vortex string. After dividing Eq.(8) by $ and 
using Eqs.(5), (6) and (7), we get 

1 f R" 1 

-i(Ni + iN 2 ) • R = - — (Ni + iN 2 ) • |-|^2 + 2V In |$| + 2i(VArg $ + eA) J 

(9) 

Identifying the real and imaginary parts separately, we get the equation of motion 
for vortex strings, 

* = <*'-kw x {*■'*<?** * +eA >) -i w x (-|^ +2Vln < 10 > 

where a is undetermined. One convenient parameterization is such that & = 
and so R • R' = 0, in which case a vanishes. This is the equation which gives 
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the velocity of a vortex string once the field configuration $ is known at a given 
moment. However, it is not a dynamical equation of motion in ususal sense because 
we do not in general know the field $ in terms of a given configuration of vortices. 
Sometimes, we can make a good guess of $ and will get some insight into the 
vortex dynamics. Furthermore, the above equation might be a good start for 
understanding two vortices crossing each other. This equation would be a sort of 
nonrelativistic limit of the vortex string equation in the relativistic field theory [3] , 
as the term proportional to R x R'/m is negligible in that limit. 

Let us study the implications of Eq.(lO) on the vortex dynamics. The field <£> 
would not change much over the scale £ near a vortex and would be determined at 
a point near a vortex by the shape of the vortex and other vortices and by other 
factors. The local fluid velocity near the vortex would be a sum of the average 
induced velocity and the current which goes around the vortex string. Naturally, 
the vortex motion would be determined by the average induced velocity, 

u = — (VArg$ + eA) (11) 
m 

In Eq.(lO), the second term in the right hand side would lead to the velocity of 
the vortex string at given point is R = u, if we neglect the last term. This is the 
standard equation for vortices in incompressible fluids. The above expression for 
the average velocity near the vortex is clearly better than the naive expression, 
(V8 + eA)/m which is very singular at the vortex position. 

The term proportional to R' x R" is interesting. The self-induction contribition 
to the fluid velocity (11) near the vortex, as we will see, has such a term with much 
larger coefficient. Thus we can regard the term proportional to R' x R" as a small 
correction to the induced velocity near the vortex. 

The last term in Eq.(lO) is relevant when the charge density has a spatial 
dependence, which is possible because the fluid is compressible. To understand its 
physics concretely, let us consider a vortex effectively moving on a two dimensional 
(x,y) plane where the charge density is forced to decrease along a direction, say 
along the +x direction. We are imagining here, for example, a three dimensional 
fluid on a container whose bottom is not leveled and so the effective two dimensional 
density decreases along the +x axis, and we assume a vortex is straight along the 
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vertical direction. Eq.(lO) tells us that it would move along the +y direction. The 
physics behind this motion is clear. The fluid rotation around the vortex would 
introduce more fluid to the — y direction from the —x direction and takes out more 
fluid from +y direction to the —x direction. Thus, the fluid pressure difference 
would push the vortex to the +y direction. However, we can expect the net fluid 
velocity is zero far way from vortex as the above phenomena is simply a continuous 
exchange of position among the vortex and the neighboring fluid. This effect can 
in principle be tested in a experiment of a rotating superfluid in a bucket whose 
bottom is tilted. This phenomena looks somewhat similar to an ascending oil drop 
in water by buoyancy. 

Since our fluid is compressible, there are sound waves due to the fluctuations, 
\& = a/po + 5pe l ° . Linearized equation for the fluctuation implies that bp = —559/g 
and the sound speed is 



^ (12) 
m 

From Eq.(lO), we see the coupling of sound waves to vortex strings is primarily due 
the fluid current fluctuations 89 in the long-wave length limit. The coupling via 
the density fluctuation term would become weaker by inverse of the wave length. 

To understand the interaction between sound wave and vortices, it is best to 
use the dual formulation of the original theory. Since the dual formulation has been 
derived many times in past [5], we will summerize briefly. We start with variables 
p, 9 with the measure [dpd9]. We introduce an auxilliary vector field, I such that 

J [dl] exp J d 4 x ^I 2 -I-W)}= p 3 / 2 exp J d 4 x £^(V#) 2 j (13) 

where irrelevant numerical factors are dropped. The phase variable is given as 
a sum of the singular part 9 and the nonsingular single- valued part rj. Because 
[d9] = [d9drf\, we can now integrate over the single valued 77, getting the charge 
conservation constraint p + V ■ I = 0. This can be solved by introducing the anti- 
symmetric tensor field B^ v such that (p, I) M = e fJjVpa d ]y B pa /2. The dual Lagrangian 
would be then written in terms of B pv and 9. 
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The coupling between vortices and would be then e^ v pa d p 9d v B pa / 2 . This 
can be written better by introducing the topologically conserved antisymmetric 
tensor current for vortices, 

= e pupa d p dj 

= 2nJ2 [ da(R»R': - RfK)5 3 (x-R a (t,a)) (14) 

a J 

where Ra = (t,~R. a (t,a)) is the position for the a-th vortex. This current is con- 
served identically, d p K pv = and leads to a conservation charge J s dT, 1 ■ K 01 on 
any surface S without boundary, which counts the net number of vortices pierc- 
ing the surface. Since the charge density is near constant outside the vortex string 
core region of the scale £, the current conservation becomes df6 = in that outside 
region. Then, we can use Eq.(14) to get the Biot-Savart law 
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= ^E/^x^ as) 



We emphasize that the solution is valid only outside the vortex core region. From 
this one can get Arg <E> near a vortex and can calculate the induced speed u. It is 
well known [9] that the above expression for a single vortex can be approximated 
by 

A R'xC A N A , l R'xR" , (Tlt . , 1R v 
* ( ' C ' a) = \R'\( B ( B 2 |Rf ( /e) ( } 

near a vortex position. Here the large length cut off is given as L and the current 
is calculated at the point C, A C A ~ £ 2 - The first term in the right-hand side would 
be the current around the vortex string. The second term is <9jArg $ and is fixed 
by the local self-induced current. The second term is in general much bigger than 
the R' x R" term in Eq.(lO). Thus Eq.(16) gives us the $ field outside the vortex 
core region. Then Eq.(lO) allows us to calculate the vortex velocity. This is how 
the cutoff £ comes in the calcuration naturally. If we have considered magnetic flux 
vortices in an extremal type II superconductor, the natural value for the cutoff L 
would be the penetration length A which is much longer than the coherence length 



7 



The field strength of B pv is H upa = d v B pcj + d p B av + d a B vp . The dual La- 
grangian becomes 

mg^ (V^i 23 ) 2 g (H ,2, i R ^ fl71 

The measure for the path integral changes from [dpd9] to [H 12 -J dB pu dR a ]. There 
will be also the gauge fixing term for the antisymmetric tensor field. (If there is a 
gauge coupling, we can integrate over the gauge field A p , getting additional terms 
in the dual Lagrangian.) The above expression is exact. Due to the second term in 
Eq.(17), the energy of a vortex does not diverge at the vortex position only if the 
charge density vanishes at the position at the vortex. Naively, there seems that a 
vortex equation would arise from the above dual Lagrangian under the variation 
of the vortex position <5R. However, the coefficient H pup of the variation vanishes 
at the vortex position and so there is no additional equation. The reason behind 
is of course that the dynamics of vortices is completely determined by the field 
dynamics around vortices. Classically, the antisymmetric tensor field is related to 
the original current by 

= l - e ^P°H vpa (18) 

In the absence of vortices, the energy of the system would be minimum when 
the expectation of #123 is the background charge density po- The small fluctuations 
around this background are sound waves and described by a simple Lagrangian 



in the long wave length limit. Here we have shifted the antisymmetric tensor field 

When we consider the small perturbations on the field configuration of moving 
vortices, we can regard the charge density #123 = p to be almost constant outside 
the vortex string core region of the transverse radius £. With the cutoff near 
vortices of length scale £, the dual Lagrangian can be written as 



Cd = _ f {Hl23 ~ po)2 " \ B » vK * v (20) 



Here we neglect the terms higher order in spatial derivatives or density fluctuations. 



8 



This is the effective Lagrangian for vortex strings and the sound wave. We can 
go further by choosing the gauge for B^ u so that d l Bij = and diB$i = 0. First, 
there will be decoupling between the Boi and Bij in the Lagrangian. We can now 
integrate over the auxiliary variables Bq{. The variation of Bqi leads to the Gauss 
law 

m 
Po 



'-d]B 0i = 2tt V / di? a <J(x - R a ) (21) 
i „ J 



We substitute Boi satisfying Eq.(21) to Eq.(20). We shift the spatial components 
of the antisymmetric tensor field as 

Bij = \ P oe ljk x k + ^e ijk C k (22) 

where C satisfies the gauge V x C = and can be regared as the fluctuation 
around the uniform charge density. 

After that we end up with the effective Lagrangian for vortices and sound waves 
with a cut off scale £ 

a a,b 

+ J rf 3 x|ic 2 -|(V-C) 2 | + 27ry^^y'rfR a -C(t,R a )xR a 

(23) 

where the gauge fixing condition V x C = is assumed. The first two terms has 
been studied before [6]. The above action is the effective Lagrangian for vortices 
interacting with sound waves. It has been known that the vortex equation has 
many interesting time dependent solutions like, Kelvin waves, Hasimoto solitons, 
and smoke rings. Perhaps, the above action will allow us to calculate the emission 
rate of sound waves by these solutions. This effective action might be also useful in 
understanding the superfluid turbulence. We can quantized the reduced effective 
action by the path integral formalism or the canonical formalism. That could allow 
us to calculate the energy levels of the excited vortex string loops and the transition 
amplitudes between the energy levels by the phonon emission. 
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Since superfluid is compressible, there could be supersonic shock waves which 
interact with vortex strings. In addition, vortex strings could move faster than the 
sound speed, still being nonrelativistic. It would be interesting to understand this 
supersonic phenomena. The dual action (17) might be useful there. 

Finally, it would be interesting to get an effective action for magnetic flux 
vortex strings in extreme type II superconductors, interacting with vector gauge 
bosons of small mass. This effective action probably approaches the effective action 
(23) in the limit where the electromagnetic coupling constant e goes to zero. 
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